Relevance of quantum fluctuations in the Anderson-Kondo model 
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We study a localized spin coupled to an Anderson impurity to model the situation found in higher 
transition metal or rare earth compounds like e.g. LaMnOs or Gd monopnictides. We find that, 
even for large quantum numbers of the localized spin, quantum fluctuations play an essential role for 
the case of ferromagnetic coupling between the spin and the impurity levels. For antiferromagnetic 
coupling, a description in terms of a classical spin is appropriate. 



I. INTRODUCTION 

Transition metal oxides show a fascinating complex 
behavior in their electronic properties i This complexity 
stems from the interplay between the formation of nar- 
row 3d-bands leading to a delocalization of these states 
on the one hand and the local part of the Coulomb in- 
teraction between the 3d-electrons tending to localize 
themi While compounds of the early 3d elements like 
e.g. LaTiOa, which typically accommodate one 3ci elec- 
tron, can at least qualitatively be described in terms of a 
one-band Hubbard mode&i^, materials involving higher 
transition metal elements like LaMnOs or TlSr2Co05 re- 
quire the use of a model including the full Zd shell. In 
particular, to understand the magnetic properties and 
the frequently occurring metal-insulator transitiongi one 
has to take into account the interplay between density- 
( "Hubbard C/") and exchange- type ("Hund's J") contri- 
butions to the local Coulomb interaction. Note that sim- 
ilar features can also be found in compounds involving 
higher rare earth elements, for example the rare earth 
pnictides. 

A particularly interesting example is Lai_2;Caa;Mn03. 
Besides is complicated phase diagram comprising a large 
variety of paramagnetic and magnetically ordered metal- 
lic and insulating phases one finds a colossal magneto- 
resistance (CMR),^ In this cubic perovskite the five- fold 
degenerate id level is split by crystal field into three-fold 
degenerate t2g, which have the lower energy, and two- fold 
degenerate Cg states. These states have to be filled with 
A—x electrons, nominally yielding a metal even for x — Q. 
However, taking into account the local Coulomb interac- 
tion, three of these electrons will occupy the t2g-states 
forming an 5 = 3/2 high-spin state due to Hund's cou- 
pling, which interacts ferromagnetically with the electron 
occupying the Cg states. Ignoring the Coulomb repulsion 
among the electrons in the eg subsystem, one encounters 
the well-known double-exchange modelSi, which has been 
extensively studied as suitable model for manganites (see 
e.g. references in Ref. In most of these investiga- 
tions, however, the t2g-spin was approximated by a clas- 
sical moment to allow the use of standard techniques like 
e.g. quantum Monte-CarlcSi2*i2iiLiaiiii4 (QMC). With- 
out such a replacement, one is typicall y re stricted to low- 
order diagrammatic techniques (Ref. [I^ and references 
therein) . 



A more realistic treatment should of course also in- 
clude the local Coulomb interaction within the Cg sub- 
system. Such a model has been proposed recently^ and 
studied in the framework of the dynamical mean-field 
theorjii^ (DMFT). Again, the t2g-spm had to be replaced 
by a classical moment to allow the solution of the DMFT 
equations with QMC. 

In this paper we want to study the validity of approx- 
imating the quantum spin by a classical object. To this 
end, we investigate the simplest possible model, viz a 
quantum impurity model consisting of a local orbital with 
interacting charge degrees of freedom coupled to a non- 
interacting host and a spin. 

While such a model surely cannot access every aspect 
of the physics of the corresponding lattice model, it is the 
basic ingredient in a DMFT calculation and thus under- 
standing its fundamental properties is of importance to 
properly interpret results obtained in a DMFT calcula- 
tion. Moreover, although such a calculation will focus on 
local dynamics only, one can obtain at least qualitative 
results about possible ordered phases, too^ To this end 
it is viable to obtain a feeling how the additional spin 
will modify local charge and spin properties. 

We employ Wilson's numerical renormalization 
groupiii^ (NRG) to solve this model. This technique 
allows to treat the model in the whole parameter regime 
and in particular identify small energy scales if present. 

The paper is organized as follows. In the next section 
we present the model and briefly review its properties for 
a classical spin in the limit of vanishing Coulomb inter- 
action. The presentation of our results follows in section 
III. The paper closes with a summary and discussion. 



II. THE MODEL 



The simplest model that allows to obtain an idea how 
the coupling to an additional local spin-degree of freedom 
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influences the properties of correlated electrons is 

a, (3 

The first three terms represent a conventional single- 
impurity Anderson modelji^iSli while the last contribu- 
tion introduces an additional spin degree of freedom 
which couples to the local states via an exchange in- 
teraction. Similar models have been studied with var- 
ious techniques in connection with double quantum 
dots.^-^i^^i^^i^*^!^^!^^ However, in these cases the additional 
impurity was represented by a correlated charge degree of 
freedom coupled via a hopping. In the limit of half-filling 
and weak interimpurity hopping, this system maps to our 
model with 5 = 1/2 and vanishing antiferromagnetic J^. 
This limit shows interesting physics on its own,^^-^^ for 
example two-stage Kondo screening or quantum phase 
transitions. However, these models neither do allow for 
ferromagnetic Jk, large and possibly anisotropic Jk nor 
spins S > 1/2, which are the particular cases we are in- 
terested in here. 

Of course our model ^ does not fully represent the 
situation found in e.g. LaMnOs since it lacks the orbital 
degrees of freedom. On the other hand, NRG calculations 
for multi-orbital models are extremely expensive^^ and 
we believe that as far as the qualitative aspects are con- 
cerned this simplification will not substantially modify 
the validity of our observations for the more complicated 
model. 

For a classical spin S, the model ^ can be solved 
exactly for U = O.^^.SQ ^j^g result for the single-particle 
Green function of the d states is 



2 \z - ed + j Ao + Jk z- ed + iAq - Jk J 

(2) 

where Aq = irNpV'^ and we assumed |5p = 1 and a flat 
conduction density of states (DOS) of inflnite width and 
value Np. The resultant DOS Pd{^) = — -^^mGdii^J + iS) 
shows two peaks of width Aq centered at Cd ± Jk ■ Note 
that one can view this result as spin-averaged DOS of 
the SIAM at [/ = in a magnetic fleld Jk and that the 
result is independent of the sign of Jk ■ 

III. RESULTS 

A. Classical limit 

Neither for a quantum spin nor for a classical spin and 
t/ > an exact solution exists. However, it is tempting 
to extend the above interpretation for a classical spin in 



the following way: Solve a standard SIAM in a magnetic 
fleld of Jk and average over the resulting spectra for 
spin up and down. Note that this procedure again leads 
to results that do not depend on the sign of Jk- 

Since the additional spin enters only on the local level, 
we use the standard NRG algorithmiSi^iiS to solve the 
impurity problem and calculate physical quantities. To 
obtain reliable spectra in a magnetic fleld, we further- 
more employ the technique proposed by Hofstetteri^^ 
Since a discretization of the energy axis introduced in 
the NRG leads to discrete spectra, a broadening must be 
introduced to obtain smooth results for dynamical quan- 
tities. Again, we follow the standard procedure herei^ 

As example we present in Fig. Q] calculations for the 
SIAM in a magnetic field, averaged over the spin direc- 




FIG. 1: Impurity DOS for SIAM in a magnetic field, with 
Ising spin 5 = 1/2 and "classical" spin S = 10. For the latter 
results for both antiferromagnetic and ferromagnetic coupling 
are shown. The inset displays an enlarged view of the region 
around uj = 0. The dashed black curve is the result for S — 
and B = 0. 

tion, together with results for the SIAM with additional 
Ising respectively classical spin. The parameters for the 
SIAM are U/ttAq = 5.3 and = -U/2, i.e. particle- 
hole symmetry. The NRG discretization parameter was 
A = 2.5 and we kept 1000 . . . 4000 states per iteration de- 
pending on the size of the local spin. The NRG-spectra 
finally were broadened with a width b = 0.6. 

The black curve in Fig. ^ was obtained from a calcula- 
tion with a magnetic field B — 8 ■ 10~^Ao, the red curve 
with a local Ising spin S* = 1/2 and coupling Jk = 4_B, 
the blue and magenta curves with local spin S' = 10 
and couplings Jk = ±2i?/10 to simulate the classical 
limit S" — > oo. The values of Jk were scaled such that 
B = JKSdS. The inset shows an enlarged view of the 
region around a; = 0; for comparison, the dashed black 
curve represents the spectrum for S — 0, B = 0. 

While the position of the Hubbard bands at higher en- 
ergies coincides for all curves, the distribution of spectral 
weight comes out different for > compared to the 
case with magnetic field and S = 0. This deviation is a 
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purely numerical effect related to the the differences in 
the distribution of spectral weight in the Hubbard bands 
for the discrete NRG spectra with and without applied 
magnetic field and the broadening introduced to obtain 
smooth spectra. 

In the region around lu = the calculation with Ising 
spin and magnetic field coincide perfectly, yielding a 
splitting of the original Kondo resonance at B = Jk = 0, 
as expected. However, the "classical limit" with S* = 10 
differs considerably. Depending on the sign of J^, either 
a Kondo resonance (ferromagnetic coupling) or a gap (an- 
tiferromagnetic case) emerges at w = in addition to the 
splitting of the original Kondo peak. Apparently, even for 
such a large value of S the effect of quantum fluctuations 
is still prominent. However, the corresponding energy 
scales are considerably reduced compared to Jk — and 
we expect the results to converge to the anticipated one 
as S' — > oo. 



B. Quantum spins: 5*= 1/2 

Let us now turn to the discussion of the effects of a local 
quantum spin on the low-energy properties. In Fig.|2]the 
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FIG. 2: Entropy and effective impurity moment for different 
values of Jk as function of T/T^^ . Full (dashed) lines denote 
antiferromagnetic (ferromagnetic) coupling. 

impurity contribution to the entropy (upper panel) and 
the effective impurity moment /i^g := T ■ Ximp (lower 
panel) for different values of Jk for a local spin S — 



1/2 are shown as function QiT/T^\ where T^' denotes 
the Kondo scale for the system with 5' = 0. The SIAM 
parameters were U — 6.47rAo at particle-hole symmetry. 
The NRG discretization, number of states kept etc. were 
chosen as before. The full lines in Fig. El represent results 
for antiferromagnetic coupling Jk i the dashed lines those 
for ferromagnetic coupling. 
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Decreasing Jk from Jk = to some antiferromagnetic 
I Jk| < results in the scenario depicted by the full red 

curves in Fig. [51 Around the temperature screening 
occurs as in the normal Kondo effect, resulting in a situa- 
tion that resembles a free spin S = 1/2 again. For a much 

lower Tk csc exp ^— a/d Jk|/Tk^'')^ , a second screen- 
ing takes place to the ground state with S — Q. This 
latter can be viewed as a conventional Kondo screening 
of the local spin by the fully formed local Fermi liquid. 
Hence, the factors appearing in the formula for Tk, 
representing the effective bandwidth respectively DOS at 
the Fermi level of the quasi-particles of the local Fermi 
liquid playing the part of the "conduction states" . This 
effect has been observed before by several authorsSSiSSiSi 
and baptized two-stage Kondo screening. When Jk be- 
comes of the order of the Kondo screening is re- 
placed by the formation of a local singlet with an energy 
scale ~ |Jk|. 

For ferromagnetic coupling Jk > 0, on the other hand, 
the proper fixed point is - as in the conventional Kondo 
model - the local-moment one with residual entropy ln(2) 
and effective moment 1/2. Again, from the dashed curves 
in Fig. 121 one can distinguish two regimes. For Jk < t'^^ 

we observe screening on the scale of , the additional 
local spin effectively behaving like a free spin all the way 
down to T = 0. However, for Jk > T^^ the impurity 
first forms a local spin triplet (entropy ln(3) and moment 
2/3). This moment is then partially screened at a reduced 
Kondo scale Tk ^ ; energy scales and physical prop- 
erties will behave as in the conventional underscreened 
Kondo modeli^^*^ 
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FIG. 3: Local DOS for different values of Jk as function of 
cj/Aq. The insets show the spectra for > in a semi- 
logarithmic plot. 

The behavior discussed previously is refiected in the 
local density of states (DOS) depicted in Fig. As al- 
ready noted in the classical limit, the coupling to the 
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additional spin leads to a corresponding shift of the up- 
per Hubbard band, which however is for larger | Jk| more 
pronounced for antiferromagnetic exchange. In this case 
on also nicely sees the two-stage screening at | Jk | < Tj^^ 

and the formation of the local singlet at | Jk | > T^'' (inset 
to upper panel of Fig. , suppressing the Kondo screen- 
ing. Here, one always finds a gap in the DOS at w = 0, 
which size is set by Tk- For ferromagnetic coupling, on 
the other hand, the inset in the lower panel of Fig. O 
proves that the screening resonance remains intact, but 
shows a width decreasing with increasing Jk- 



with increasing Si. Note that even for comparatively 
large 5*/ =3/2 the influence of quantum fluctuations is 
still pronounced and appears in an energy regime that 
may still be of experimental relevance. 

The differences become even more pronounced if we in- 
crease the local Coulomb repulsion to U/Aq ~ 10, which 
lies in the intermediate-coupling regime with a Kondo 

scale Tj^^V^o ~ 0.07. Since thus Jk > T^'' we do not 
expect two-stage screening here. The results for oth- 
erwise same model parameters are collected in Fig. |S1 
Note that due to the choice of Jk the Hubbard bands do 



C. Quantum spins: General S 

How does the behavior discussed in the previous sec- 
tion change with increasing spin quantum number Sj or 
more precisely, do we recover a "classical" result for large 
enough Sjl Let us start with a discussion of the weak- 
coupling results shown in Fig. ^ The calculations were 
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FIG. 4: Local DOS for fixed ^Si[Si + 1)|Jk|/Ao = 1/3 as 
function of w/Ao in the weak-coupling regime U/Aq — 1 for 
different values of local spin Si. The upper panel collects the 
results for antiferromagnetic coupling, the lower panel those 
for ferromagnetic. The insets show the spectra for w > in a 
semi-logarithmic plot. 

done in the weak-coupling regime U/ Aq — 1 for fixed 
value \/Si{Si + 1)|Jk|/Ao = 1/3 to achieve the same 
classical energy scale for all Si. We did calculations up 
to Si — 10, which, according to our results in section 
nil Al we expect to be already very close to the classical 
limit. Indeed, for Si = 10 (green curves in Fig.^J we do 
find almost identical behavior for Jk < and Jk > ex- 
cept for extremely low temperatures. Note however, that 
even for this large value for Si the DOS for uj/Aq < 0.1 
for Jk < does not reach the full unitary limit due to 
quantum fluctuations. 

The differences are more dramatic for small values of 
Si. As expected, for Jk < there occurs a Kondo screen- 
ing with an energy-scale T}<i{Si) decreasing exponentially 




FIG. 5: Local DOS for fixed y^Si{Si + 1)|Jk|/Ao = 1/3 
as function of u/Ao in the intermediate-coupling regime 
17/Ao — 10 for different values of local spin Si. The up- 
per panel collects the results for antiferromagnetic coupling, 
the lower panel those for ferromagnetic. The insets show the 
spectra for > in a semi- logarithmic plot. 

not move with increasing Si. Furthermore, in all cases 
Si > additional structures appear at u ~ ±| Jk|. The 
low-energy behavior, however, is markedly different for 
Jk < (upper panel in Fig. [SJ and Jk > (lower panel 
in Fig. El). Below uj/Aq < 0.07 « Tj^\ the former case 
always develops a (pseudo-) gap due to the formation of 
a singlet between the local degrees of freedom, while the 
latter tends to recover a Kondo-resonance for the total 
spin, again with exponentially decreasing Tk{Si). Ow- 
ing to the universal behavior of the conventional Ander- 
son model for U / ttAq > 1 ( "Kondo regime" ) we actually 
expect the behavior observed here to be generic in this 
parameter regime. Quite obviously, while in the weak 
coupling regime the "classical limit" is reached already 
for moderate values of Si, one has to be very careful when 
dealing with the strongly correlated regime. 



IV. SUMMARY AND CONCLUSION 

In this paper we presented calculations for an extended 
Anderson impurity model, where the local charge degrees 
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of freedom in addition couple to a localized spin. The 
motivation to study such a model is based on the obser- 
vation that in a variety of transition metal or rare earth 
compounds the complex local orbital structure can be 
split into a localized spin, which can take large values 
S ^ 1/2, coupled via Hund's exc;hange to a more delo- 
calized set of possibly also correlated states. A particular 
example sm-ely is the famous LaMnOs. 

The solution of this model for different regimes of 
model parameters was accomplished by using Wilson's 
NRG, which provides accurate and reliable results for 
thermodynamics and dynamics and is able to resolve ar- 
bitrarily small energy-scales that may appear in the prob- 
lem. The findings can be summarized as follows: Even for 
comparatively large localized spin Si = 10, we still ob- 
serve the influence of quantum fluctuations on the prop- 
erties of the impurity charge degrees of freedom. These 
effects become more pronounced when these charge de- 
grees of freedom are correlated themselves as to be ex- 
pected for example in LaMnOa. Depending on the ratio 
\Jk\IT^\ where is the Kondo temperature for the 
model without additional spin, different regimes can be 
identified, which in contrast to the classical prediction do 
markedly depend on the sign of Jk- 

Thus, for the solution of correlated lattice models with 
such an additional spin degree of freedom within DMFT 
one has to be likely careful when using the approxima- 



tion of a classical spin, even when Si is comparatively 
large. Moreover, the expected physics can be read off 
our results right away, at least for half-filling. Due to 
the reduced Kondo scale for Hund's type or ferrmomag- 
netic coupling, we expect a corresponding reduction of a 
critical U for a Mott-Hubbard transition. On the other 
hand, for antiferromagnetic exchange coupling the corre- 
sponding ?7c = at T = 0, because the forming of a local 
singlet immediately leads to an insulating state. 

Quite interesting are also the magnetic properties of 
the system. Again, we may anticipate from the impurity 
calculations that antiferromagnetism is still the prevail- 
ing magnetic order, but in the vicinity of certain com- 
mensurate band fillings we can also expect ferromagnetic 
order from a corresponding RKKY exchange. These in- 
verstigations are currently in progress. 
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